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We calculate the spectrum of a beam of atoms output 
from a single mode atomic cavity. The output coupling uses 
an internal state change to an untrapped state. We present 
an analytical solution for the output energy spectrum from a 
broadband coupler of this type. An example of such an output 
coupler, which we discuss in detail uses a Raman transition 
to produce a non-trapped state. 
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As a result of recent experiments in which a Bose Ein- 
stein Condensate (BEC) has been produced in the lab 
there has been considerable interest in coupling the 
atoms in a BEC out of a trap. This could produce a con- 
tinuous, coherent, directional beam of atoms - an atom 
laser beam f^dTl"!. While initial experiments have suc- 
ceeded in coupling atoms out of a BEC by changing the 
internal state of the atoms to a non-trapped state |l2],[l3| , 
there is still much to be understood about the output 
beam. In this paper we present an analytical solution for 
the output energy spectrum of atoms in a single trapped 
mode coupled to free space by a change of internal state. 
Our analysis is based on the atom field input-output the- 
ory presented by Hope juj. We discuss the dependence 
of the spectrum on output coupling strength, and relate 
these findings to the MIT atom laser experiment |l2],[l3| . 

In a BEC a large number of bosonic atoms are cooled 
into a single energy eigenstate of a trap. This is an im- 
portant step towards producing a monoenergetic beam 
of atoms. Nevertheless we still have the problem of how 
to coherently couple the atoms out of such a trap in a 
way that preserves their monoenergetic nature. 

There are many ways in which atoms can be coupled 
out of a trap. The simplest method is to turn off the 
trap Q. The result of rapidly turning off the trap is to 
reproduce the BEC wavefunction in free space. In par- 
ticular, the wavefunction momentum width is conserved. 
As a result, the atoms have the corresponding range of 
energies in free space and the monoenergetic nature of 
the original BEC is lost. Fortunately, energy conserving 
output coupling is possible. One example is quantum 
mechanical tunneling of atoms through the trap walls. 
This is the atomic analogue to the use of partially trans- 
parent mirrors on an optical laser. Such a process has 
been considered in a model of an atom laser proposed by 
Wiseman ||. It would be difficult in practice, however, 
to use tunneling to produce sufficient fluxes of atoms due 
to the exponential dependence of the tunneling rate on 



the trap potential barrier. 

Another approach to the output coupling problem 
would be to change the internal state of the trapped 
atoms to an untrapped state. Experimentally such a 
method has been used by implementing radio-frequency 
pulses to induce spin flips on trapped atoms in a BEC 
p^ , ^3| . Furthermore the use of Raman transitions as 
a method of output coupling has been suggested [ p"l"| |. 
Raman transitions have a number of advantages. A Ra- 
man transition can have an extremely narrow linewidth 
so that lasers can be tuned so as to only couple atoms 
from a particular trap mode, due to energy conserva- 
tion. Moreover when Raman beams are oriented so that 
they are counter propagating, they provide a momentum 
kick of size 2hk. This could be used to provide direc- 
tionality to the atomic output beam if atoms were sup- 
ported against gravity, for instance in a hollow optical 
fiber |§-0. 

We model here an output coupler based on change 
of state, focusing initially on the specific case of a Ra- 
man output coupler which uses two lasers tuned to a 
two-photon resonance to couple atoms between an initial 
atomic state, and a final atomic state. There is a third, 
excited, atomic state which mediates the Raman transi- 
tion. We assume that each of the lasers is far detuned 
from single photon resonance. In this far detuned limit 
we can adiabatically eliminate the third state to produce 
an effective two level Hamiltonian. In this Hamiltonian 
we ignore the energies of higher atomic modes of the trap. 
Initially these other modes are empty as we assume all 
the atoms are condensed in the ground mode. Ignoring 
these higher energy modes for later time is valid for very 
narrow linewidth Raman lasers which are only on res- 
onance with the ground trap mode. This ensures that 
higher modes do not become populated by atoms in the 
output state transferring back into the initial state at 
later times. In addition population of other modes is 
suppressed by Bose enhancement of transitions into the 
ground mode |ll]]. We also ignore the effects of atom- 
atom interactions. The resulting effective Hamiltonian is 
then of the form 
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Here, the single trap mode is described by the creation 
operator, and is coupled by the Raman lasers to a 
continuous spectrum of external modes described by cre- 
ation operators, bl. hu>x (H0J2) is the energy of the trap 
(output) atomic state. Tlloq is the ground state trap en- 
ergy, m is the mass of the trapped atoms. hk\L and 
hk2L are the momenta of the two lasers inducing the 
Raman transition, with frequencies uj\l and u>2L respec- 
tively. Thus %{k\L + &2l) is the total momentum kick re- 
ceived by atoms making the Raman transition, fii (CI2) is 
the Rabi frequency of the transition between the trapped 
(output) state and the excited state which mediates the 
Raman transition. Ai and A2 are the detunings of the 
two Raman lasers from the excited state. We have as- 
sumed these are large in adiabatically eliminating the up- 
per level. If the lasers are tuned close to the two-photon 
resonance, Ai w A2. ip(k) is the momentum space wave- 
function of the ground mode of the trap. T is a coupling 
strength, given here in terms of the Rabi frequencies and 
single photon detuning. 

The form of the Hamiltonian, Eqs. (|| - f|), is valid 
in the more general case of an arbitrary output cou- 
pling through state change involving a single mode sys- 
tem coupled to a continuous spectrum of external modes. 
In this more general case Kuq gives the energy of the 
trapped atoms, and hui k gives the energy of the free 
atoms. The coupling strength is more generally defined 
through n(k,t) = T2K,'(k,t) where n'(k,t) describes only 
the shape of the coupling and is normalised to unity. 
The form of n'(k, t) for a general interaction describing a 
change of state is n(k, t) = ip*(k — k ) Q . Here, ip(k) is 
the ground state momentum space wavefunction of the 
single mode system and ko describes a possible fixed mo- 
mentum kick applied to the atoms in the state change 
process. In the following we discuss the Raman coupling 
case, given by Eqs. (|| - ||) for definiteness. The results, 
however, are valid for a general output coupler in the 
regime where the coupling strength, T, and the energies 
huiQ and Tiujk are suitably defined. 

We are interested in the output energy spectrum, 
(blbk) which is the mean population density of the contin- 
uum of free space momentum eigenstate modes, labelled 
by the momentum hk. We obtain this by solving the 
Heisenberg equations of motion for the operators, b k (t). 
In general, such a solution is difficult to obtain, however 
recently Hope jl4j has presented a solution in terms of 
inverse Laplace transforms. Using these solutions, the 
output spectrum, in the case where initially the external 



modes are empty is given by 
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The final equality holds for the case when the lasers are 
tuned to the two photon resonance in free space, which 
we assume here. L and C~ l are the Laplace transform 
and inverse Laplace transform respectively. 

We present an analytic solution for the spectrum in the 
limit of broadband coupling. For simplicity, we consider 
the case where the total momentum kick from the Raman 
lasers is very small. That is we assume k\L ~ —k2L- This 
is analogous to the MIT output coupling experiments in 
which the atoms receive a negligible momentum kick in 
changing state Jl2|,[l3| . We also assume that the coupling- 
function n(k, t) is broad. The shape of n(k, t) is given by 
the ground state momentum wavefunction of the trap, 
ijj(k). We consider here a harmonic trap, with a gaussian 
ground state of standard deviation a k in wavenumber 
space. We can calculate an exact value for £(/')(s) from 
the definition given in Eq. ([ll]) , however we must simplify 
£(/')(s) in order to evaluate Eq. (|l0|). In the regime 
where Im(s) << ha\jra we can approximate £(/')(s) by 
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Using this approximation to calculate M k it) is equivalent 
to discarding high (> ha 2 /m) frequency information in 
the Laplace transform space. As we increase the width 
of our coupling in momentum space, given by a k , our 
solution for M k (t) becomes valid for increasingly high 
frequencies. For an infinitely broad coupling our expres- 
sion becomes exact, and is equivalent to the form of the 
general broadband coupling discussed by Hope |Q. Us- 
ing the above expression for C(f')(s) we find the inverse 
Laplace transform, Mfc(t) to be 
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where we have defined 0Jk — hk 2 / (2m) and A& = u>k — ujq. 
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The function L L 2 (a;) is a Laguerre polynomial, Erf is the 

2 

error function and a, (3 and 7 are the roots of the equation 
s 3 + iujQS + Tc^/i = 0. 

Fig. [I] shows the behaviour of |Mfc(t)| 2 as a function of 
Uk and time after we turn on the output coupling inter- 
action. Initially \Mk(t)\ 2 is small, and for short enough 
times, arbitrarily broad in k-space. Initially |Mfc(i)| 2 
agrees with the perturbative solutions presented by Hope 
jnj. For longer times, we can see that the spectrum 
reaches a stable shape. For very large values of the cou- 
pling strength, the long time limit becomes very broad in 
k-space. As a result, the shape of the output spectrum, as 
given by Eq. (|^) , simply reflects the momentum distribu- 
tion of the cavity wave- function, tp(k). As a result there 
is no narrowing of linewidth in momentum space. The 
recent MIT experiments fl2l , [l3| are an example of an out- 
put coupling with an extremely large coupling strength. 
In these experiments a short, 5/is RF pulse was used to 
couple atoms out of a BEC, making a pulsed atom laser. 

We consider here a continuous coupler, turned on at 
time t = 0, and examine the resulting long time spectrum 
in the external modes described by b\. We observe in 
Fig. [IJthat for longer times |-Mfc(i)| 2 narrows into a sine 
function centered about the trap ground state frequency, 
ujq. Eventually |Mfc(i)| 2 reaches a stationary state with a 
lorentzian like profile as shown in Fig. [j]. This longtime 
behaviour is given by 
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where 7 is the particular solution to the cubic discussed 
above, given by the expression 
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The longtime expression for Mu (t) , Eq. ( jig ) contains two 
terms. The first of these terms dominates in the case of 
small T, while the second dominates for very large T. 
As a result, the long time spectrum has two distinct be- 
haviours depending on the strength of the coupling. We 
consider the case of slow coupling (small T) initially. In 
this case, the long time expression for Mk(t) is dominated 
by the first term in Eq. (|16|) above, and the resulting long 
time spectrum is given by 
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A plot of the long time spectrum, Eq. (|18| ) as a function 
of u!k is presented in Fig. ||for various coupling strengths. 
Fig. H shows that for increasing coupling strength the 
linewidth of the long time spectrum increases. The val- 
ues for r chosen correspond approximately to values of 
Raman laser Rabi frequencies, f^i ps 2n x 50 kHz and 
Q, 2 ps 2tt x 1.6 MHz and detuning, Ax ps 2vr x 2.5 GHz 
similar to values presented in However, much 

smaller or larger coupling strengths can be achieved by 
suitably adjusting the intensities of the lasers and their 
detunings. The figures assume a trap with ground state 
frequency luq = 2n x 123 s _1 , typical of magnetic traps 
for ultra cold atoms |is|] . A ground state gaussian with 
width <jfc ps 10 m has been assumed, which corre- 
sponds to a position space wavefunction of size of the 
order of 2/im. This value of Ofc corresponds to a width in 
u>k space of a LUk ps 10 4 s -1 . 

For each of the graphs shown in Fig. ||, the lorentzian 
like spectrum is centred about u>o, the ground state fre- 
quency of the single mode trap, with the width of the 
spectrum dependent on the strength of the coupling as 
mentioned above. In all cases, however, the linewidth is 
much less than that which would be obtained if the trap 
was rapidly turned off, that is <j Uk ps 10 4 s _1 . We see from 
Eq. (|l8|) that the distribution isn't exactly lorentzian due 
to the presence of lju in the second part of the denomi- 
nator. However for large ujq the spectrum is well approx- 
imated by a lorentzian distribution of width IFcI/y/oJo. 

We have already noted that for large coupling rates, 
the width of the longtime limit of |Affc| 2 , and hence of 
the longtime spectrum is increased. When T is very 
large, \Tc\/-y/uJo >> a^ kl the width of Mfc(i) becomes 
large compared with n(k,t) and the spectrum becomes 
dominated by the cavity momentum spread ip(k). As 
a result, for sufficiently fast coupling (large T) the out- 
put spectrum changes significantly from the lorentzian 
shape considered above, and instead reflects the momen- 
tum spread of the cavity. This is shown in Fig. ||. For 
very large T the spectrum is centred about zero, and falls 
away exponentially in Uk space, as required for a gaussian 
distribution in momentum space given by ip(k). 

We have shown that the longtime spectrum from an 
output coupler based on state change depends on the 
strength of the output coupling. For very strong cou- 
pling, the output spectrum is given by the cavity spec- 
trum, and is very broad in momentum space. The spec- 
trum is then centered about the zero of momentum when 
there is no net momentum kick from the lasers. As 
the strength of the coupling is reduced, however, the 
long time linewidth is correspondingly reduced. For 
small coupling strengths the final linewidth is effectively 
lorentzian, centred about the energy of the cavity with a 
linewidth proportional to the coupling strength T. 
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FIG. 1. Plot of |Mfc(£)| 2 as a function of cjk and time for 
t — Os to t = 5s, and LJk ranging from 762s _1 to 783s _1 about 
the single mode trap frequency, luq ~ 772s -1 . T = 1.8x 10 3 s -2 



FIG. 2. Plot of the long time behaviour of (b k bk) as a func- 
tion of oj fc for various coupling strengths, T = 10 4 s~ 2 (dotted 
line), T = 3 x 10 4 s~ 2 (solid line) and T = 5 x 10V 2 (dashed 
line) . 



FIG. 3. Plot of the steady state behaviour of {b k b k ) as a 
function of uj k for the large coupling limit (r « 10 13 s~ 2 ). 
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